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Abstract. We investigate the dynamics of bound states of two interacting 
particles, either bosons or fermions, performing a continuous-time quantum walk 
on a one-dimensional lattice. We consider the situation where the distance 
between both particles has a hard bound, and the richer situation where the 
particles are bound by a smooth confining potential. The main emphasis is on 
the velocity characterizing the ballistic spreading of these bound states, and on the 
structure of the asymptotic distribution profile of their center-of-mass coordinate. 
The latter profile generically exhibits many internal fronts. 
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1. Introduction 

Quantum walks have witnessed an upsurge of interest in parallel with the developments 
of quantum algorithms and quantum information (see [1, 2, 3] for reviews). Two 
different types of quantum-mechanical analogues of classical random walks have been 
investigated. Discrete-time quantum walkers [4, 5, 6, 7] possess, besides their spatial 
position, a finite-dimensional internal degree of freedom, referred to as a quantum coin. 
Both spatial and internal degrees of freedom jointly undergo a unitary dynamics. 
Continuous-time quantum walkers [8, 9, 10] have no internal degree of freedom. 
Their dynamics is governed by some hopping operator on the underlying structure. 
Despite these differences, continuous-time quantum walks can be viewed as a limit 
of discrete-time quantum walks [11]. Both types of quantum walks exhibit many 
similar properties. Their main characteristic feature is a fast ballistic spreading. The 
typical distance traveled by a quantum walker grows linearly in time, as opposed to 
the diffusive spreading of a classical random walker. 

If two or more quantum walkers are simultaneously present, the combined effects 
of interactions, quantum statistics and entanglement give rise to novel features which 
have no classical counterpart. The Anderson localization of two interacting quantum 
particles in a random potential has attracted much attention [12, 13, 14]. More 
recently, dynamical features of the quantum walks performed by two or more entangled 
or interacting particles have been the subject of numerous theoretical studies (see 
e.g. [15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]). Several experimental groups have also 
studied the quantum walk of entangled pairs of magnons [26] and of photons in various 
integrated photonics devices [27, 28, 29, 30]. 

Here we investigate the continuous-time quantum walk of bound states of two 
identical bosonic or fermionic particles. The main emphasis is on the distribution 
profile in the center-of-mass coordinate of the bound states, including the dependence 
of the ballistic spreading velocity on model parameters, and the generic presence of 
many internal ballistic fronts. The setup of this paper is as follows. In section 2 we 
revisit in a pedagogical way the continuous-time quantum walk of a single particle. We 
analyze the distribution of the position of the particle, emphasizing its dependence 
on the initial quantum state. We also show that allowing the particle to hop to 
second neighbors may yield a novel effect, namely the appearance of internal ballistic 
fronts in the position distribution, besides the usual extremal fronts. We then address 
similar questions in the more challenging situation of the quantum walk performed by 
bound states of two identical particles. The same formalism allows one to deal with 
bosonic and fermionic bound states, as they are respectively described by even and odd 
functions of the relative coordinate between both particles. We consider bound states 
generated either by imposing a hard bound on the distance between both particles 
(section 3) or by a smooth confining potential (section 4). In both situations the 
emphasis is on the distribution of the center-of-mass coordinate and on the presence 
of internal ballistic fronts. Section 5 contains a discussion of our findings. 

2. Quantum walk of a single particle 

This section serves as a self-contained pedagogical introduction to the main concepts 
emphasized in the rest of the paper. We consider the continuous-time quantum walk 
performed by a single particle on a one-dimensional lattice. Most of the features 
underlined in this section have been studied, or at least mentioned, in several earlier 
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works in the context of discrete-time quantum walks [20, 31, 32, 33, 34, 35, 36, 37]. One 
of the purposes of this section is to demonstrate that the analysis of the continuous¬ 
time quantum walk problem is much simpler, and can therefore be worked out in a 
more systematic fashion. 

More specifically, we first revisit the simple quantum walk where the particle hops 
to nearest neighbors only, focusing our attention onto the asymptotic distribution of 
the walker, including its dependence on the initial state. We then turn to a generalized 
quantum walk, where the particle hops to further neighbors as well. We show that 
allowing hopping to second neighbors gives rise to internal fronts in the distribution 
profile. Hopping to second neighbors is known to have drastic physical consequences 
in many other situations. One celebrated example is graphene (see [38, 39]), where 
hopping to second neighbors breaks the chiral symmetry between both sublattices. 

2.1. Simple quantum walk 

The framework of the simple quantum walk is the usual one of the tight-binding 
approximation, where the particle hops to nearest neighbors only. Throughout this 
paper, we use dimensionless units. The wavefunction ipnit) = (nj'i/'(t)) of the particle 
at site n at time t obeys 

i = ^"-1-1 (^) -I- '^ri-i{t). (2.1) 

The dispersion law between wavevector (momentum) q and frequency (energy) to and 
the corresponding group velocity v therefore read 

(jj{q) = 2cosq, v{q) = uj'{q) = —2smq, (2.2) 

where the prime denotes a derivative. 

Initial state localized at the origin 

Suppose that at time t = 0 the particle is localized on a single site, taken as the origin 
of the lattice: ■0ra(O) = ^n,o- In Fourier space, ^{q,0) = 1, hence f:{q,t) = 
and so 

Mt) = r p- = i-” J„(2t), (2.3) 

where the Jn are the Bessel functions. 

Figure 1 shows the probabilities against position n at time 

t = 50. These probabilities exhibit various regimes of behavior at large n and t. 
They take appreciable values in the allowed region which spreads ballistically with the 
maximal velocity 

V =max la;'(( 7 )l = 2. (2.4) 

Furthermore, they exhibit sharp fronts near n = ±2t, and they decay exponentially 
in the forbidden region beyond these fronts. These asymptotic results are classics of 
the theory of Bessel functions [40, 41], whose physical meaning has been underlined 
in [9]. These results can be readily obtained the saddle-point method, which will be 
used later in other situations. 
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n 


Figure 1. Probability profile \tpn{t)\‘^ at time t = 50 for the simple 
quantum walk with a particle initially located at the origin. Vertical blue 
lines: theoretical positions of the ballistic fronts. Black curve: asymptotic 
distribution (2.7) set to scale. 


• Allowed region (|n| < 2t). In this region, the oscillatory behavior of the wave- 
function at long times can be derived as follows. Setting 

n = -2tsing* (|g*|< 7 r/ 2 ), (2.5) 

the integral in (2.3) is dominated by two equivalent saddle points at 9 = g* and 
q = TT — q^, (modulo 27r). We thus recover the well-known asymptotic form of 
Bessel functions: 

V Trt cos g* 

By averaging out the oscillations in the probabilities |^„(t)p, one arrives [9] at a 
smooth distribution f(v) for the ratio v = n/t in the long-time limit, i.e., 

f{v) = (k| < 2). (2.7) 

The above result can be alternatively obtained by folding the uniform distribution 
on the Brillouin zone [—tt, tt] by the dispersion curve of the group velocity: 
V = —2 sing. 

• Transition region (|n| Ri 2t). The distribution (2.7) becomes singular as the 
endpoints of the allowed region are approached (n —>■ ±2<). The vicinity of these 
endpoints corresponds to the transition region in the theory of Bessel functions. 
Setting \n\ =2t-\- we have 

J„(2t) = i”V^„(t) r: t“^/^Ai(z), (2.8) 

where Ai( 2 :) is the Airy function. The probabilities |V'n(i)P therefore exhibit sharp 
ballistic fronts [9], whose height and width respectively scale as t~'^/^ and t^/^. 

• Forbidden region (|n| > 2t). In this region, the exponential fall-off of the 
wavefunction at long times can be derived by evaluating again the integral in (2.3) 
by the saddle-point method. Setting |n| = 2f cosh0 with 0 > 0, we obtain 

^ — 2t{0 cosh sinh 0) 


J„(2t) = 'Flfnit) 


'JAnt sinh 9 


(2.9) 
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Arbitrary initial state 

For an arbitrary initial state i/'n(0)) we have 


Jo 


Whenever the initial state is reasonably localized, in the sense that tpn{Q) decays 
fast enough with distance |n|, the time-dependent wavefunction exhibits qualitatively 
the regimes of behavior described above in the allowed region, ballistic fronts, and 
forbidden region. 

Let us analyze more carefully the allowed region (|n| < 2t). With the 
definition (2.5), the integral in (2.10) is now dominated by two inequivalent saddle 
points at q = qi, and g = tt — g* (modulo 27r). We thus arrive at 


/Itt cos g* 


9 *, 0 ) e' 


i(ng* — 2t cos g* +7r/4) 


-h $iTr - g*, 0) eH"(’^-9*)+2‘cos9,-V4)^ _ (2.11) 


Averaging out the oscillations in the above expression yields the following formula 
for the locally coarse-grained probabilities: 


\i’n{t)\ 


2 _ l^(g*: 0 )P + 1^(74- g,^, 0 )P 

47 rt cos g* 


The limit distribution f(v) of the ratio v = n/t is obtained by folding the distribution 
on the Brillouin zone with density |^/’(g)p/(27r) by the dispersion curve v = —2sing. 
This line of thought has been used in [35, 36], and more thoroughly in [37], for discrete¬ 
time walks. 

For definiteness, let us consider the case where the initial wavefunction is spread 
over three consecutive sites, i.e., 

'*/’n(0) = a<5n,l + bdnfi cSn,-l (joj^ -|- Jl)]^ -|- Jcj^ = 1). (2.13) 


We have 


f:{q,0)=ae bce" 


(2.14) 


fjn{t) = i "(6J„(2t) -t i(aJ„_i(2t) - cJ„+i(2t))). 

The asymptotic result (2.12) reads explicitly 

...2 1-I-4A cos^ g* - 2B sin g* 

IV'n(t)| « --, 

Zirt cos g* 

with the definition (2.5), and where): 

A = Re(ac), i? = Im(6(a — c)). 

The limit distribution f{v) of the ratio v = n/t reads therefore 

^ I + A{2 - + Bv ,, , 

f{v) = -Sj==- (bl<2). 

TTV 4 — 


f{v) = 


(2.15) 


(2.18) 


I The bar denotes complex conjugation. 
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Analogous expressions have been derived in [5, 6, 35, 36, 37] for various models of 
discrete-time quantum walks. The distribution (2.18) generically has inverse-square- 
root singularities at the endpoints of the allowed region {v = ±2), corresponding to the 
ballistic fronts. The initial state only affects the numerator (see (2.7)). This lack of 
universality, namely the ever-lasting memory of the initial state, is a genuine quantum 
feature (it is absent for the classical walker with sufficiently localized initial state). 

The hrst moments of the distribution (2.18), 

{v) = 2B, (u2) = 2(1-A), (2.19) 

match the asymptotic growth at long times of the exact expressions 

= ( 2 . 20 ) 

of the position moments [9], namely 

(n(t)) = \a\'^ - \c\‘^2Bt, 

{n^(t)) = \a\^ + |c|^ -f 2 Im(6(a -f c))t + 2(1 - A)t^. (2.21) 

The inverse-square-root singularities at the endpoints of the allowed region 
(v = ±2) are generic but not entirely universal. There are indeed special initial 
states such that either one or both endpoint singularities are absent from the limit 
distribution (2.18). Figure 2 shows the probability profiles |^/>„(t)p at time t = 50 in 
the following two special cases of initial states. 

1. For a = l/-\/2, b = and c = 0, hence A = 0 and B = 1/2, we have 

n 

and so the left front is absent (see figure 2, left). 

2. For a = c = l/-\/2 and 6 = 0, hence A = 1/2 and B = 0, we have 

jl-n 

M) = -^(J„_i(2t) - j„+i(2t)), (2.24) 

f{v) = (kl < 2), (2.25) 

and so both fronts are absent (see figure 2, right). 


{Jn-l{2't) Jn{2tf), 


( 2 . 22 ) 


This kind of non-generic behavior of the distribution prohle has also been 
observed for discrete-time quantum walks, both with the usual two-dimensional 
internal state [32, 20] and with a more exotic three-dimensional quantum coin which 
may lead to a localization phenomenon, in the sense that part of the probability stays 
forever at a hnite distance from the particle’s starting point [33, 37]. Reference [37] 
also contains analytical predictions for /(u), somehow similar to (2.23) and (2.25), in 
many special situations. 
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n n 

Figure 2. Probability profiles \'tpnit)\'^ at time t = 50 for the simple 
quantum walk with the two exceptional initial states described in the text. 
Vertical blue lines: theoretical positions of the ballistic fronts. Black curves: 
asymptotic distributions (2.23), (2.25) set to scale. 


2.2. Generalized quantum walk 

Novel phenomena occur in generalized quantum walks, where the particle may hop to 
further neighbors. Let us consider the minimal generalization where hops are limited 
to first (nearest) and second (next-nearest) neighbors, with respective amplitudes 1 
and g. The wavefunction of the particle at site n at time t now obeys 

i ^ = 'f’n-i-lit) + 'f’n-lit) + g (V'n+2(0 + '^n-2{t)) ■ (2.26) 

at 

We have therefore 

oj{q) = 2{cosqgcos2q), 

uj'{q) = v{q) =—2{smq2gsin2q), 

uj"{q)= — 2{cosq4:gcos2q). (2.27) 

The above dispersion law ^(< 7 ) is invariant under the transformation g —> —( 7 , q —> 9 + 7 r, 
u! —^ —oj. We may therefore restrict the analysis to the domain g > 0, without any 
loss of generality. 

Suppose the particle is launched from the origin: ipniO) = Sn,o- Thus ip{q,Q) = 1, 
and so 

V;„(t)= / ^ Qdnq-u{q)t) ^ ^2.28) 

Jo 27r 

Some observables have a smooth dependence on the amplitude g. This is 
especially the case for the position moments (see ( 2 . 20 )), which read 

(n 2 (^))= 2 (l+V)^^ 

(n^(t)) = 2(1 -b 16g2)t2 -b 6(1 -b l&g"^ -b 16/)^^ (2.29) 

and so on (odd moments vanish identically by symmetry). 

The presence of hopping to second neighbors however introduces a novel 
qualitative feature. Let us for the time being adopt a general standpoint, and consider 
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the wavefunction (2.28) in the regime of long times. Evaluating the integral by the 
saddle-point method, as we did in section 2 . 1 , we obtain 


g 


Qi(nq-U)(q)t) 

A/27ila7^(^’ 


(2.30) 


where the sum runs over the solutions q of the saddle-point equation 
n 


t 


= v = Lu'{q). 


(2.31) 


By averaging out the oscillations in the probabilities again predict that v 

has a smooth distribution 




1 


27r ^ |w"(g)| 


(2.32) 


in the long-time limit. In full generality, the above distribution is again obtained by 
folding the uniform distribution on the Brillouin zone by the dispersion map of the 
group velocity v = w>'{q), in the sense that dg/(27r) = f(v)dv holds formally. 

The distribution (2.32) has generically inverse-square-root singularities at all 
extremal values of v, in correspondence with wavevectors q so that dv/dq = uj"{q) 
vanishes. It is always singular at the endpoints of the allowed region {v = ±E), as 
the maximal velocity V is necessarily an extremal value. The distribution (2.32) may 
however also have internal singularities within the allowed region, which were absent 
in the simple quantum walk considered in section 2 . 1 . 

Internal ballistic fronts of that kind have been observed in two generalizations of 
the usual discrete-time quantum walk [31, 34]. Reference [31] investigates a quantum 
walk subjected to M independent quantum coins acting cyclically. If M is large, the 
dynamics exhibits a crossover between classical random walk at short times {t ^ M) 
and quantum walk at long times (t^M). In the latter regime the distribution of the 
quantum particle exhibits an array of equally spaced ballistic peaks, whose number 
grows as M/2, the distance between any two consecutive peaks being Au = ^/2/M. 
The model studied in [34] is closer to ours in its spirit. It consists of a discrete-time 
quantum walk where hops up to distance j are allowed, whose dynamics is rigidly 
dictated by (2j -|- l)-dimensional Wigner rotation matrices. Here too, the distribution 
profile exhibits 2 j -|- 1 equally spaced ballistic peaks. 

In order to pursue the analysis, let us specialize to the continuous-time quantum 
walk with hops to first and second neighbors only (see (2.26), (2.27)). This minimal 
example is already too complex to allow one to turn (2.32) to an explicit expression. 
The location of the ballistic peaks can nevertheless be predicted as follows. The second 
derivative uj"{q) vanishes for 


cosq± = 


-I± (I -h 12852)^/2 

16^ 


The corresponding values V± of v are given by 


(2.33) 


2 -1-^32052+ 2048^4 ±(1-^12852)3/2 

“ 12852 


(2.34) 


The allowed region always spreads ballistically with the maximal velocity V+. The 
smaller velocity V- may also play a role, depending on the strength of 5 : 
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• If the amplitude g is small enough {g < g^ = 1/4), the situation is qualitatively 
similar to that of the simple quantum walk, studied in section 2.1. We have 
indeed cosg_ < —1 and 0 < cosq+ < 1, so that only matters, and so f{v) is 
only singular at the endpoints ±V+ of the allowed region. 

• If the amplitude g is large enough {g > gc = 1/4:), both and q- matter. As a 
consequence, the distribution f{v) has two singularities at the endpoints ±V+ and 
two internal singularities at the smaller values ±1A. The wavefunction exhibits 
four ballistic fronts: two extremal ones, propagating at the maximal velocity 
(±V+), and two internal ones, propagating at a smaller velocity (±VT). 

• In the borderline case {g = gc = 1/4), we have cosq- = —1, i.e., <?_ = tt and 
V- = 0. The dispersion curve exhibits an unusual quartic behavior near q = tt: 

o 4 

w(g)«-2 + ^ (e = (7-TT0). (2.35) 

This anomalous dispersion right at g = gc has two consequences. First, the 
distribution f(v) has a singularity at u = 0, of the form 


f{v) 


1 

67r|up/3 


(2.36) 


This central singularity with exponent —2/3 is stronger than the generic ones, 
whose exponent is —1/2. Second, the wavefunction at the origin exhibits an 
unusually slow fall-off: 


Mt) 


V2r(5/4) eh3‘/2-V8) 
nt^/^ 


(2.37) 


This t decay is slower than the generic t decay exhibited e.g. by the Bessel 
function Jo(2t) (see (2.3)). 


Figure 3 shows the dependence of the front velocities against g lor g > 0. The 
larger velocities ± 1 + of the extremal fronts (blue curves) describe the endpoints of 
the allowed region for all g. The smaller velocities ±VT of the internal fronts (red 
curves) exist for g > gc = 1/4 only. At g = gc we have V+ = 3-\/3/2, while V- takes 
off as {32'/6/9){g — gc)^^^■ At large g, both velocities grow linearly in g with the same 
slope, as V± « 4(/ ± This non-trivial pattern of front velocities is richer than the 
periodic arrays of equally spaced peaks observed in [31, 34]. 

Figure 4 shows the probabilities | 7 /„(t)p at time t = 50 for a particle launched 
at the origin and two values of g. For g = 1/2 (left), the probability profile exhibits 
four fronts. For g = gc = 1/4 (right), the probability profile exhibits three fronts. The 
central one at the origin corresponds to the anomalous singularity (2.36). 


3. Two-body bound states: hard bound on distance 

We now turn to our main subject: the quantum walk performed by a bound state of 
two identical particles propagating coherently along a one-dimensional lattice. The 
main emphasis will be on the asymptotic distribution of the center-of-mass coordinate 
of the bound state, on the velocity characterizing its ballistic spreading and on the 
structure of the distribution profile, which generically exhibits many internal fronts. 
To the best of our knowledge, these matters have only been addressed so far in two 
papers [21, 25]. In order not to interrupt the lengthy developments of sections 3 and 4, 
we postpone the discussion of those earlier works to section 5. 
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Figure 3. Front velocities of the quantum walk with hopping to first and 
second neighbors (see (2.26)), against the amplitude g. External (blue) 
curves: velocities ±1/+ of the extremal fronts (endpoints of allowed region). 
Internal (red) curves: velocities ±V1 of the internal fronts for g > Qc = 1/4. 



n n 


Figure 4. Probabilities for the quantum walk (2.26) at time 

t = 50. Vertical blue lines: theoretical positions of the ballistic fronts. 
Left {g = 1/2): four fronts at ±V+ and ±l/_. Right {g = gc = 1/4): 
two fronts at ±V+ and a central one at the origin, corresponding to the 
anomalous power-law singularity (2.36). 


Here again, considering continuous-time walks will allow for a more thorough and 
systematic investigation of the problem. Bound states obtained by imposing a hard 
bound £ on the distance between both particles are dealt with in this section, whereas 
those generated by a smooth confining potential will be considered in section 4. The 
same formalism will allow one to deal with bosonic and fermionic bound states, as 
they are respectively described by even and odd functions of the relative coordinate m 
between both particles. 
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We denote by ni = n + m and n 2 = n the abscissas of two identical particles on the 
lattice, where to = ni — n 2 is the relative coordinate. We impose a hard bound £ on 
the distance |to| between both particles, and so to is restricted to the 2£ + 1 values 
TO = —£, ...,£. Figure 5 shows the positions of the particles in the (ni,n 2 ) plane for 
£ = 2. Full symbols denote the allowed configurations of the particles. Links between 
symbols show the allowed hops of any of the particles to a neighboring site. 




¥ 4.0 ¥j.i 


Vn -2 V„Ll 


• i 


Figure 5. Plot of the positions of the particles in the (ni,n 2 ) plane for 
£ = 2. Full symbols denote the allowed configurations of the particles. 
Links between symbols show the allowed hops of any of the particles to a 
neighboring site. 


The time-dependent wavefunction 

i’n,m{t) = ((ni,n2)|f/'(0) = {{n -£ m,n)\'if{f)) 
obeys the equation 

. d'0n,m(O 


dt 




“f —l,m+l (0 V’n.m-t-l (f)5 

with boundary conditions V'n,±(^+i) = 0- 


(3.1) 


(3.2) 


3.2. Bosonic and fermionic spectra 
A basis of plane-wave solutions to (3.2) reads 

V;„,„(t) = edn.mq+mp-u.{p,q)t) ^ (3 3^ 

where the momenta p and q are respectively conjugate to the relative coordinate 
TO = ni — n 2 and to the center-of-mass coordinate 


ni -b n2 TO 

ncm =-^- =n+—. 


(3.4) 


The resulting dispersion relation has a product form [25]: 
q) = dcospcos 


(3.5) 
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• Bosonic states are described by even functions under the exchange of rii and n 2 ', 
they are obtained by adding the plane waves (3.3) for p and —p: 

(t) = eiK„,9-^(p.9)t) cos mp. (3.6) 

The relative momentum p is quantized by the condition cos(^ + \)p = 0. It 
therefore takes the t' + 1 values 


(B) 

Pk = 


(^ + g)' 

i+l 


{k = 0, 


(3.7) 


and so 


(<?) = 4 cos 


{k + |)7r 

e + 1 



(3.8) 


The bosonic dispersion curve thus consists of £ +1 branches, with group velocities 




—2 cos 


(fc + 4)7r 
£+1 



(3.9) 


• Fermionic states are described by odd functions of m, obtained by subtracting 
the plane waves (3.3) for p and —p: 

sinmp. (3.10) 

As fermionic particles cannot cross each other in one dimension, the range of 
the fermionic wavefunctions (3.10) can be restricted to the sector ni > n 2 , 
i.e., m = !,...,£. The relative momentum p is quantized by the condition 
sin(^ + l)p = 0. It therefore takes the £ values 


and so 


(F) 

Pk 


kn 

J+l 


(fc = 1,.. .,£), 


/p\ 

^k (?) = 4 cos 


kiT 

7Tl 



(3.11) 


(3.12) 


The fermionic dispersion curve thus consists of £ branches, with group velocities 


JF) 




(3.13) 


Figure 6 shows the bosonic and fermionic spectra for a maximal distance £ = Q. 

/p'l /"P'1 

The 7 bosonic frequencies uj\. {q) (black) and the 6 fermionic frequencies uj). {q) (red) 
are plotted against q/r:. 


3.3. Ballistic spreading 

Let us now investigate the asymptotic behavior of the wavefunction ipn.mit), starting 
from an arbitrary initial state, where both particles are located in the vicinity of the 
origin. The detailed analysis of the one-body problem performed in section 2 allows 
one to draw the following picture. 

The various components 'f’n.mft) of the wavefunction spread ballistically, i.e., their 
extension in the center-of-mass coordinate ncm ~ n grows asymptotically linearly in 
time t and symmetrically with respect to the origin. The ipn.mit) generically exhibit 
sharp fronts near n = Vkt, where the front velocities 14 are the stationary values of the 
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0 0.2 0.4 0.6 0.8 1 


q/7t 


Figure 6. Plot of the bosonic and fermionic spectra for a maximal distance 
£ — 6, against g/rr. Black: bosonic frequencies Red: fermionic 

frequencies {Q)- Blue: band boundaries (a; = ±4cos(g/2)) (see (3.5)). 


group velocity, corresponding to the boundaries of the Brillouin zone. Setting q = zLtt 
fn (3^), (3.13), we get 


(fc+ i)7 


kvr 


= 2 cos ■ 


1 


(k = 0,...,£), 
(k = l,. 


(3.14) 


In particular, all the components of the wavefunction take appreciable values only 
in the allowed zone defined by |n| < Vt, where the maximal velocities 


1/(^1 = 2 cos ■ 


= 2cos- 


2 (^+ 1 )’ ’ ^ 

respectively correspond to setting k = 0 and fc = 1 in (3.14). Both maximal velocities 
approach the free value 14 = 2 in the limit of a large confining size (£ ^ 1), with two 
different correction amplitudes, i.e.. 


a; 2 - 


= 


C(B) 


a; 2 - 


C7(f) 

~£^'‘ 




(3.16) 

(3.17) 


For a generic initial state localized in the vicinity of the origin, the various 
components tpn,m(t) of the wavefunction will exhibit, besides two extremal ballistic 
fronts at n = or n = £ — \ internal fronts in the bosonic case (for 

£>2) and £ — 2 internal fronts in the fermionic case (for > 3). 


S.f. Continuum limit and corrections 

Our bound-state problem owes its non-triviality and its richness to the fact that 
particles live on a lattice. In the continuum limit, the dynamics of the center-of-mass 
coordinate and of the relative coordinate are exactly decoupled, as a consequence of 
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Galilean invariance. As soon as the line is discretized into a lattice, this Galilean 
invariance is broken. This phenomenon has already been underlined for discrete-time 
interacting quantum walkers [21]. Its consequences in the context of diluted Fermi 
gases and of the BCS-BEC crossover have also been discussed recently [42]. 

In the present situation, too, the dynamics of the center-of-mass coordinate ncm 
and of the relative coordinate m are expected to decouple as the continuum limit is 
approached, where Galilean invariance should be restored. It is worth investigating in 
a quantitative way how this decoupling takes place, starting from the lattice dispersion 
curves (3.8) and (3.12). For that purpose, we introduce the lattice spacing a, and define 
the following dimensionful quantities, denoted by capital letters. The half-width of the 
potential well is conveniently defined as L = (£ -{- l)a, while the conjugate momentum 
to the center-of-mass coordinate is Q = q/a, and finally the continuum energy reads 
E = 

To leading order as a —>■ 0 (i.e., in the vicinity of the center of the Brillouin zone), 
both dispersion relations (3.8) and (3.12) yield 

Ek^^+Sk, (3.18) 

where the Sk are the energy levels of a free particle in a potential well of width 2L in 
the appropriate sectors, i.e., 



The expression (3.18) conforms to what we expect from the continuum theory: the 
total energy Ek of the compound system is the sum of the kinetic energy Q^/d of the 
center-of-mass motion and of the energy Sk of a bound state in the relative coordinate. 

The corrections to the leading-order result (3.18) can be derived by recasting (3.8) 
and (3.12) in terms of the dimensionful quantities introduced above, and expanding 
in powers of a. The hrst correction thus obtained, 

Ek = Sk — -I- 24Q‘^ek Ibe^) H-, (3.20) 

already demonstrates that the coupling of both degrees of freedom by the lattice 
structure affects the energy spectrum in a non-trivial way. 


3.5. A case study: fermionic state with maximal distance I = A 


To close this section, let us investigate in detail the dynamics of a fermionic bound 
state with maximal distance t = A. This is the smallest £ for which generic dynamical 
behavior is observed.§ The front velocities (3.14) read 


vP=-vP 
yf) = - yf) 


2 

2 


TT -x/b -I-1 

27r - 1 

COS — = T = --- 


(3.21) 


where r is the golden mean and r its reciprocal. 

Suppose the two fermions are launched from sites 0 and 1 at time t = 0. Then 
V'ra,m(0) = 6n,odm,i and the initial mean value of the center-of-mass coordinate is 


For £ = 3 , there exists a central anomalous front with zero speed. 
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(’T^cm(O)) = (V'(O)l^cmlV'(O)) = 1/2- We expect there will be two extremal ballistic 
fronts near n = ±Tt, and two internal ones near n = ±Tt. 

The most efficient way of solving the differential equations (3.1) with prescribed 
initial values consists in performing a spatial Fourier transform (with q being conjugate 
to n) and a temporal Laplace transform (with s being conjugate to t). The Fourier- 
Laplace transforms thus defined obey the equations 

isf^iiq,s) = (e"‘« + 1 )^( 2 (<?, s) + i, 
isV' 2 (g, s) = (e‘« + 1)^1 (g, s) + (e"‘« + l)V^3(g, s), 
isipaiq, s) = (e‘« + l)V^ 2 (g, s) + (e"“* + 1 )V^ 4 (<?, s), 
isipiiq, s) = (e‘« + 1 )V^ 3 ( 9 , s)- (3.22) 

This linear system can be readily solved. Introducing the notation 

7 = 2cos|, (3.23) 

we obtain 


■*Ai(< 7, s) 
i’3{q,s) 


s(s^ + 27 ^) 
A(g,s) 

(e‘« + 1)^5 
A{q, s) 



with 


A{q, s) = {s — i7r)(s — i7T)(s + i7T)(s + i7r). 

The inverse Laplace transform of (3.24) can be taken first. We have e.g. 


V’i(g, s) = — -j= ( --1--1-H-) ) 

2v5 \ s ~ 17 T s + 17 T s — 17 T s + 17 T / 


and so 


Mq,t) = — ’ 


The inverse Fourier transform can then be taken using (see (2.3)) 

f ^ gi(n9+a.cos(g/2)) ^ (_ 1 )"J 2 „(^), 

Jo 27r 

r" ^ gi((n+l/2)g+xcos(9/2)) ^ i(_ 1 )" . 

Jo 27 r 


We thus obtain 


V’„,i(t) = —^(TJ2„(2Tt) + rJ2„(2rt)), 

v5 

'^n,2{t) = i--(^2n+l(2Tt) + J2ri+1 ( 2 rt)) , 

'f’n,3{t) = -p-(J2n+2(2rt) — J2n+2(2Tt)), 

V5 

f—Ij" 

tpnA(J) = i ^ ^ (rJ2n+3(2rf) - rJ2„+3(2rt)). 

V5 


(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 


Figure 7 shows the four probability profiles |V’n.m(i)P thus obtained at time t = 50. 
Abscissas are shifted from n to n + (m — l)/2 = ncm — (f^cm(O)), in such a way that 
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the plotted profiles are exactly symmetric. The four profiles exhibit the same global 
features, although they differ in their detailed structure. As predicted, they exhibit the 
same ballistic fronts, two external ones at n « ±rt and two internal ones at n ~ ±rt 
(see (3.21)). The probabilities are larger within the internal fronts, and considerably 
smaller in the wings of the allowed zone, i.e., between the internal and the external 
fronts. This phenomenon is quite generic; it could already be observed in the left 
panel of figure 4. 



-100 -50 


0 

n 


50 100 



n+1/2 




n+1 


-100 -50 0 

n+3/2 


50 100 


Figure 7. Probability profiles (m = 1 to 4) of the fermionic 

bound state with maximal distance £ = 4 at time t = 50, with initial state 
t^n,m(0) = (5n,0(5m,i, agaiust n-l- (m — l)/2 = ricm — (ncm(O)). Vertical blue 
lines: theoretical positions of the ballistic fronts. 


Let us now turn to the internal structure of the fermionic bound state. The 
distance between both particles takes the value m = 1,..., 4 with probability 


Pm{t) = ^ |'0n,m(i)P- 

n 

These probabilities can be evaluated from (3.29) using identities 

n 

r ^ ^ _ Mx-y)-Mx + y) 

y , J2n-i-l{x)J2n+l{y) — -^-, 


(3.30) 


(3.31) 
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which can be derived from the integral representations (3.28). Quadratic identities 
of this kind have their roots in the connection between special functions and 
representation theory [43]. We thus obtain 

Pi(t) = ^(3 + ?Vo(4rt) +t2Jo( 4?<) + 2Jo(2t) + 2Jo(2y5t)), 

P2{t) = ^(2 - JoiiTt) - Jo(4ft) + 2Jo(2t) - 2Jo{2V5t)), 

Ps{t) = ^(2 + Jo(4Tt) + Jo(4ft) - 2Jo(2t) - 2Jo{2V5t)), 

P4{t) = ^(3 - rVo(4rt) - r^MiTt) - 2Jo(2t) + 2Jo(2^51)).(3.32) 


These probabilities sum up to unity, as should be. They are even functions of t whose 
power series have rational coefficients. At short times we have 

51^ 35t® 
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P,(t) = l-2t^ + ^- 
3 It® 

P^it) = 2f - 4^4 + ^ + 


Psit) 

P4{t) 


3t^ 5t® 

'2 Y 

51® 

_l- 


In the long-time regime, these probabilities reach the stationary values 


(3.33) 


Pl=P4 


10 ’ 



(3.34) 


These numbers agree with the result (A. 14) derived in Appendix A. In the present 
situation, it is indeed legitimate to study the stationary properties of the internal state 
per se, without referring to the ballistic dynamics of the center-of-mass coordinate, 
because the compound system has a basis of factorized eigenstates (3.3). This feature 
is by no means general. 

Figure 8 shows a plot of the probabilities Pm{t) against time t. The maximal 
time 1 = 50 corresponds to the profiles shown in figure 7. The stationary values (3.34) 
(arrows) are reached after a complex pattern of damped oscillations. The envelope of 
this transient oscillatory behavior falls off very slowly as while the oscillations 

themselves follow a quasiperiodic pattern, characterized by the two incommensurate 
frequencies 2t and 2r. All the frequencies entering (3.32) are indeed integer linear 
combinations of the latter frequencies. In a generic situation, there will be as many 
incommensurate frequencies as there are positive ballistic velocities Vf or Vf 
{see (3.14)). 

The mean internal size of the bound state, 

4 

D{t) = Y, ^Pm{t), (3.35) 

m—1 

can be evaluated from (3.32) to read 


D{t) = ^(25 - (2 - 3T)Jo(4rt) -{2 + 3r) Jo(4?t) - 8Jo(2t)). (3.36) 

This quantity is plotted in figure 9. It starts from D{0) = 1 and reaches the stationary 
value D = 5/2, again after a complex pattern of quasiperiodic oscillations dying off 
as t-P\ 
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Figure 8. Probabilities Pm{t) describing the internal structure of the 
fermionic bound state (see text) against time t. Arrows: stationary 
values (3.34). 



Figure 9. Plot of the mean internal size D{t) of the fermionic bound state 
against time t. Arrow: stationary value (D = 5/2). 


4. Two-body bound states: smooth confining potential 

4-1. Generalities 

We now turn to the quantum walk performed by a bosonic or fermionic bound state 
of two identical particles generated by a smooth confining potential Wm ■ The latter is 
assumed to be an even function of the relative coordinate m between both particles, 
such that Wm —t +oo at large distances (|m| —>■ oo). With the notations of section 3, 
the time-dependent wavefunction ipn,m(t) obeys the equation 

■ d'tjjn^rnit) ^ (0 + 1 (^) +V'n+l,m-lW 

at 

H- V^n,m+l(0' (^■^) 

Let us look for a basis of solutions to (4.1) of the form (see (3.3)) 

(4.2) 
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where the momentum q is conjugate to the center-of-mass coordinate (3.4). The 


internal wavefunction obeys 

T fm+l) (4-3) 

with the same dispersive (i.e., g-dependent) hopping amplitude 

7 = 2 cos ^ (4.4) 

as before (see (3.23)). We shall also use the shorthand notation 

= (-1)”^</)™, (4.5) 

for the staggered wavefunction, which obeys 

WXm = -7(Xm-l + Xm+l) + WmXm, (4.6) 

i.e., (4.3) with the sign of 7 reversed. 


The above formalism applies to an arbitrary confining potential. Bosonic and 
fermionic bound states are described by wavefunctions (fm which are respectively 
even and odd functions of to, such that (fm —>■ 0 at large distances. Such bound-state 
wavefunctions only exist for discrete sequences of dispersive frequencies LO^\q) and 
Hereafter we investigate in detail the case of potentials growing as a power 
of distance, either linearly (section 4.2), quadratically (section 4.3), or as an arbitrary 
power (section 4.4). 

4 . 2 . Linear confining potential 

Our first example is the linear confining potential 

W^=g\m\, (4.7) 

where the amplitude g is a positive constant. The regime of most physical interest 
corresponds to small g, so that bound states have a large size and potentially a rich 
internal structure. 

Figure 10 shows the bosonic and fermionic spectra for g = 0.4. These spectra 
exhibit two very distinct regions. Within the band delimited by the blue curves 
(w = ±27 = ±4cos(g/2)), bosonic and fermionic branches are well separated from 
each other. They alternate and are strongly dispersive. Above the band, the spectrum 
consists of an infinite array of approximately equally spaced and non-dispersive 
branches. These branches appear as red, as bosonic and fermionic branches are 
superimposed to a very high accuracy. 

In order to analyze the above observations, we use (4.6) for positive values of the 
relative coordinate to, i.e., 

WXm = -7(Xm-l + Xm+l) + gmXm (w > 0). (4.8) 

Let US first forget about the constraint to > 0 and extend (4.8) to all values 
of TO. We are thus facing a tight-binding equation for a charged particle in a uniform 
electric field, the amplitude g giving the strength of the field in reduced units. The 
corresponding spectrum is a Wannier-Stark ladder [44, 45]: it consists of the infinite 
sequence of equally spaced frequencies, 

Wfc = gk, 


(4.9) 
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0 0.2 0.4 0.6 0.8 1 

q/% 

Figure 10 . Plot of the bosonic and fermionic spectra for the linear 
confining potential with g = 0.4, against q/n. Black: bosonic frequencies 
Red: fermionic frequencies Blue: band boundaries 

(a; = ±27 = ±4cos((?/2)). 


where k runs over the integers. The Wannier-Stark levels are non-dispersive, 
as (4.9) holds irrespective of the hopping amplitude 7 . The normalized eigenfunction 
corresponding to ujk reads 

Xm = Jm-k{‘2^x/g), (4.10) 

where the Jm are the Bessel functions. This eigenfunction is strongly localized around 
m = k. We have in particular, for large positive k, 

This estimate shows that the Wannier-Stark eigenfunctions hardly feel the boundary 
condition on (4.8) at m = 0, and therefore remain essentially unperturbed, as soon 
as the frequency uJk = gk exceeds a few times the bandwidth 7 . In other words, 
only finitely many lowest Wannier-Stark states, in a number of the order of 7 /g, are 
affected by the boundary condition which distinguishes between fermions and bosons. 
This explains the main observations made on figure 10. 

A more quantitative analysis of the problem goes as follows. The general solution 
to (4.8) falling off as m —>■ -boo reads 

Xm = Jm-ujjgif^llg) (^ > “1), (4.12) 

where the J^iz) are the Bessel functions with index v, which obey the identities 

z 

J,_i(2)-J,+i(z) =2J'(0), (4.13) 

where the prime denotes a derivative with respect to z. The bosonic spectrum is given 
by xi = X-i, while for the fermionic spectrum xo — 0 - We thus arrive at the exact 
quantization formulas 


( 27 / 5 ) = 0 , 


J-u>(r)/g{2xlg) = 0 . 


( 4 . 14 ) 
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The leading correction to the Wannier-Stark spectrum (4.9) at large positive k 
can be readily derived from (4.14). Skipping details, we obtain a symmetric splitting 
of the form 

^ g{k-k Ek), ui^p ^ g{k - Ek), (4.15) 


with 


Ek 


ii/gr 

fc!(fc- 1)! 




(4.16) 


where xo is the Wannier-Stark wavefunction at the origin (see (4.11)). 

Our next goal is to investigate the velocities characterizing the ballistic spreading 
of a bosonic or fermionic wavefunction in the center-of-mass coordinate Ucm, for a 
generic initial state. These velocities are given by 

y(i') =inax \u;f\q)\, =max |a;f ^(g)|, (4.17) 

k,q k^q 


where primes denote derivatives, and maxima are taken over all branches of the 
spectrum and over all q. It is clear from figure 10 that these maxima are reached 
for the lowest branch of each spectrum. Furthermore, if the amplitude g is small, 
so that the bound states have a large internal size, we anticipate that the maximal 
velocities will be close to the free value 2 and correspond to momenta close to the zone 
boundary (g —>■ tt). 

To substantiate these expectations we note that whenever g is small, typical 
wavefunctions vary slowly with m, so one can employ a continuum description. Setting 

w = —27 + d = —4 cos ^ (4-18) 

and expanding in (4.8) differences in terms of derivatives, we obtain 


Then, setting 


d^y 

7 -^ « {gm-6)x- 

dm 

x = X{gm- 6 ), A = (3^7)"^/^, 


Equation (4.19) becomes the Airy equation 
d^X 




= XX, 


whose solution decaying as x —>■ +00 is x(x) = Ai(x), the Airy function. 


(4.19) 

(4.20) 

(4.21) 


• Bosonic states obey xi = X-ii hence dy/dm = 0 for to = 0, and so xq = —\6 
obeys Ai'(xo) = 0. The lowest branch of the bosonic spectrum therefore reads 

(5(B) « ,71(52^)1/3, (4.22) 

where gi = 1.018792 ... is the opposite of the first zero of Ai'(x). 

Setting q = TT — E, the lowest branch can be expanded as 

3 

w = —4 sin ^ + (5 = —2e + ^ + • • • + A (4.23) 

Using the expression (4.22) of S, we can show that the maximal group velocity 
reads 


U(B) « 2 - C(B) 5I/2, 


( 4 . 24 ) 
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= 



= 0.551985..., 


(4.25) 


and is reached for 

£« (4.26) 

• Fermionic states obey yo = 0, and so xq = —XS obeys Ai(a:o) = 0. The lowest 
branch of the fermionic spectrum therefore reads 

(4.27) 

where = 2.338107... is the opposite of the first zero of Ai(a;). The maximal 
group velocity now reads 

« 2 - (4.28) 


C(F) 



1.029227..., 


(4.29) 


and it is reached for 

£f«(C(^))l/2 5l/4. 


(4.30) 


The above scaling results valid in the 5 <C 1 regime can be alternatively derived 
by analyzing the exact quantization formulas (4.14) in the transition region (see (2.8)). 
We have in particular 


C(F) 

C(B) 


ff \ 3/4 

= 1.864592... 


(4.31) 


Figure 11 shows plots of the maximal velocities and against 5^^^. The 
scaling results (4.24) and (4.28) at small g (straight lines) provide a good description 
of the velocities for moderate values of g. 



Figure 11. Plot of the bosonic and fermionic maximal velocities 
(black) and pT) (red), characterizing ballistic spreading in a linear 
confining potential, against Straight lines: scaling results (4.24) 

and (4.28) at small g. 
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4 . 3 . Quadratic confining potential 


Our second example is the quadratic confining potential 

W'm = gm\ (4.32) 

where g is again a positive constant. Equation (4.6) reads 

^Xm = -7(Xm-i + Xm+i) + gm^Xm- (4.33) 


The Fourier transform x{p) therefore obeys 



— (w + 27 COsp)x. 


(4.34) 


The latter differential equation is known as the Mathieu equation [46]. The body 
of knowledge on the latter equation is however of little use for the present purpose. 
Hereafter we therefore use general techniques which could be applied to any confining 
potential. 

Figure 12 shows the bosonic and fermionic spectra for g = 0.1. These spectra 
are qualitatively similar to those shown in figure 10. The existence of two very 
distinct regions, with strongly dispersive modes within the band and weakly dispersive 
branches above the band, is indeed a common feature of all confining potentials. 



0 0.2 0.4 0.6 0.8 1 

q/7I 


Figure 12. Plot of the bosonic and fermionic spectra for the quadratic 
confining potential with g = 0.1, against q/TT. Black: bosonic frequencies 
^?\q)- Red: fermionic frequencies Blue: band boundaries 

(w = ±27 = ±4cos(g/2)). 


Let us first consider the weakly dispersive branches above the band. Right at 
7 = 0 (this corresponds to g = tt, i.e., to the right end of figure 12), the eigenstates 
are strictly localized at specific sites: we have Xm = ^m,fc and ujk = Wk = gk^- If 7 is 
small and/or k is large, the wavefunction Xk±i at neighboring sites is proportional to 
the ratio 7/fc. We thus obtain the more refined estimate 


LUk ~ gk^ + 


27' 


(4fc2 - l)g' 


(4.35) 


High branches are therefore weakly dispersive, as their bandwidth scales as l/(gfc^). 
With a linear confining potential, the Wannier-Stark states were not dispersive at all. 
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In the whole region above the band, the splitting between even (bosonic) and 
odd (fermionic) frequencies can hardly be observed. This splitting can be estimated 
as follows. First, it is expected on general grounds to scale as Xo (4-16)). 
Furthermore, xo is exponentially small whenever 7 is small and/or k is large. In 
this regime, (4.33) indeed simplifies to g^k"^ — mf)xm ~ —7Xm+i for 0 < m < /c. We 
thus obtain the estimate 


1X01 


2 (7/g)" 

(2fc)! 



(4.36) 


The velocities I/i®! and characterizing the ballistic spreading of an initially 
localized wavefunction again correspond to the lowest branches of the spectra. In the 
regime of most interest where g is small, these velocities can be calculated by means 
of a continuum description. With the notation (4.18), we obtain 


7t-4 ~ (4.37) 

dm 

Setting 

m = {ghf/^x, 5 = (4.38) 


Equation (4.37) becomes the usual Schrodinger equation for a harmonic oscillator. 


d^X 2 




(4.39) 


whose eigenvalues are Ty = 2v \ (^ = 0,1,...). 

For the bosonic spectrum, the lowest branch corresponds to Sq = 1, and so 


sm „ 


(4.40) 


Using (4.23), we obtain 


]/(B) 

C(B) 


2 - 
5 

214/5 


= 0.717936... 


(4.41) 

(4.42) 


For the fermionic spectrum, the lowest branch corresponds to £i =3, and so 


« 3(57)'/", (4.43) 

« 2 - (4.44) 


C(F) 


214/5 


1.728952... 


(4.45) 


The ratio of the amplitudes is 
(7(F) 

= 3^/® = 2.408224 ... (4.46) 

Figure 13 shows plots of the maximal velocities and against 5^/®. Here 
again, the scaling results (4.41) and (4.44) at small g (straight lines) provide a good 
description of the velocities for moderate values of g. 
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Figure 13. Plot of the bosonic and fermionic maximal velocities 
(black) and (red), characterizing ballistic spreading in a qnadratic 

confining potential, against Straight lines: scaling results (4.41) 

and (4.44) at small g. 


4-4-. Confining potential with an arbitrary exponent 


We now consider the case of a confining potential of the form 

^^^=51^1“, (4.47) 

with an arbitrary exponent a > 0. 

Let us focus on the maximal velocities and characterizing the ballistic 
spreading of a localized initial state. In the most interesting regime when g is small, 
these velocities can be calculated by means of a continuum description. With the 
notation (4.18), we obtain 

7^-4 ~ (ff|mr - (5)x. (4.48) 

dm 

Setting 

m = (g/7)i/(“+2)x, S = (4.49) 

Equation (4.48) becomes the continuous Schrodinger equation 

-^ + krx = (4.50) 

whose eigenvalues (z/ = 0,1 ,...) are not known analytically in general. 

For the bosonic spectrum, the lowest branch corresponds to £q, and so 

ft! (52^“)1/(“+2 )£:o. (4.51) 

Using again (4.23), we obtain 

« 2 - 52/(a+3)^ ^4 52) 


C(B) 


Q: “t“ 3 / Ol£{) 


(Q+2)/(a+3) 


22 /(«+ 3 ) 


(a+ 2)2 


(4.53) 
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For the fermionic spectrum, the lowest branch corresponds to £i, and so 




(4.54) 

yiP) 

«2-C'(F) g2/(a-e3)^ 

(4.55) 

CiP) 

a-b3 / a£i \(“+2)/(«+3) 

“ 22/(“+3) l^(a-H2)2 j 

(4.56) 


The ratio of the amplitudes is 


C(F) 

ClB) 


(a+2)/(a+3) 


(4.57) 


The scaling laws (4.52) and (4.55) generalize (4.24), (4.28), (4.41), (4.44). They 
can also be put in perspective with (3.16), which holds in the situation where a hard 
bound £ is imposed on the distance \m\ between both particles. Let us introduce the 
length £ as the distance where the potential (4.47) equals unity, i.e., 

£ = g-^/°^. (4.58) 


In terms of this length parameter, (4.52) and (4.55) take the form 


y(B) 


CCB) 

£2a/(a+3) ’ 


yiF) 


C(F) 

£2a/(a+3) ’ 


(4.59) 


These expressions smoothly match their counterparts (3.16) in the a —>■ oo limit, 
where the potential gets infinitely steep. The exponent slowly converges to the limit 
value 2, while the amplitudes (4.53), (4.56) also converge to the limits (3.17). 

The ratio can be viewed as a universal amplitude ratio. This 

dimensionless quantity depends only on the growth exponent a of the confining 
potential (see (4.57)), increasing monotonically from 1 in the singular a —>■ 0 limit 
to 4 in the a —>■ oo (i.e., hard-bound) limit. Its values for a = 1 and a = 2 have been 
given in (4.31) and (4.46). Figure 14 shows a plot of this ratio against a/{a -b 1). 

Finally, the amplitude ratio is always larger than unity. More generally, 

bosonic bound states always have a larger maximal spreading velocity than fermionic 
bound states in the same confining potential. We shall come back to this very common 
property in section 5. 


5. Discussion 

In this work we have investigated the continuous-time quantum walk performed along 
a one-dimensional lattice by bound states of two interacting particles. The main focus 
has been on the profile of the wavefunctions describing these bound states in the 
center-of-mass coordinate, and especially on the velocity characterizing their ballistic 
spreading and on the structure of the whole profile, which generically exhibits many 
internal fronts. 

We have first revisited the problem of a single quantum walker in a self-contained 
pedagogical fashion. For the simple quantum walk where the particle hops to nearest 
neighbors only, we have concentrated onto the dependence of the distribution of the 
particle position on the initial state. This distribution profile has generically two 
ballistic fronts. Either one or even both fronts may be absent for carefully chosen 
special initial states, as this was already the case for some examples of discrete-time 
walks [20, 32, 33, 37]. For the generalized quantum walk, where the particle hops to 
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Figure 14. Plot of the universal amplitude ratio against 

a/{a + 1), where a is the growth exponent of the conhning potential. 
Symbols: results (4.31) and (4.46) for a = 1 and a = 2. Curve: data 
points obtained by means of a numerical calculation of the eigenvalues £o 
and £\ of (4.50). 


first and second neighbors with respective amplitudes 1 and g, we have emphasized 
the possible occurrence of two internal fronts in the distribution profile, propagating 
at velocities ±V-, besides the two usual external ones, propagating at velocities ±V+, 
which mark the endpoints of the allowed region beyond which the wavefunction falls off 
exponentially. The non-trivial dependence of the velocities V± on the amplitude g is in 
contrast with the equally spaced internal peaks arising for the discrete-time quantum 
walks considered in earlier works [31, 34]. 

The quantum walk of bound states of two bosonic or fermionic particles has then 
been investigated in two situations: either by imposing a hard bound i on the distance 
between both particles, or by generating the bound states by a smooth confining 
potential growing as a power of the distance. In both situations, we have focused on 
the structure of the distribution of the center-of-mass coordinate. We have investigated 
in detail the maximal velocities and characterizing the ballistic spreading 
of bosonic and fermionic bound states, as well as the many internal fronts of their 
distribution profiles. In the case of a hard distance bound (section 3), the maximal 
velocities have the simple expressions (3.15). The distribution profile generically 
exhibits i—\ (resp. £—2) internal fronts in the bosonic (resp. fermionic) case, besides 
the two extremal ballistic fronts. In the situation of a smooth confining potential 
(section 4), we have investigated in detail the cases of a linear and of a quadratic 
(i.e., discrete harmonic) potential. For all potentials of the form Wm = g\m\°‘, growing 
as a power of the distance between both particles, the maximal velocities exhibit 
scaling laws of the form (4.52), (4.55) in the regime of a weak potential {g <C 1). The 
associated amplitude ratio is universal, in the sense that it depends only 

on the growth exponent a of the confining potential. This ratio is larger than unity, 
in accord with the more general property that bosonic bound states have a larger 
maximal spreading velocity than their fermionic counterparts in the same binding 
potential. 
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Let us now compare our findings with those of two recent papers [21, 25] which 
are also devoted to the quantum walk of bound states. The latter works deal with 
full two-body Hamiltonians, along the lines of earlier investigations of the Anderson 
localization of two interacting particles [12, 13, 14]. In such models bound states 
therefore coexist with a full two-particle continuum. The situations studied in the 
present work (hard distance bound or conhning potential) only possess bound states, 
so that the investigation of the latter is made easier. As internal ballistic fronts are 
not discussed in [21, 25], we shall henceforth concentrate on the maximal spreading 
velocity V. Reference [21] deals with the antisymmetric (i.e., fermionic) sector of 
a discrete-time model, with a local interaction described by the action of a special 
coin operator whenever both particles sit at the same site. The associated coupling 
constant g is an angle. The analysis of the bound-state spectrum allows one to derive 
the maximal velocity V. This quantity (whose expression is not explicitly given in [21]) 
decreases continuously from l/-\/2 to 1/3 as g is increased from 0 to its maximal 
value of TT. Reference [25] describes the continuous-time dynamics of two identical 
particles (either bosons or fermions) generated by a Hubbard-like Hamiltonian with 
nearest-neighbor interactions. The maximal velocities are derived by means of a 
perturbative approach in the regime of very strongly attractive interactions. Both 
spreading velocities are found to fall off to zero as the inverse of the interaction 
strength, and to obey the simple relation 1/1®) = 3R*-^/ 

Even though very different models and regimes have been considered, all the 
findings recalled above are consistent with the following universal characteristics of 
the ballistic spreading of two-body bound states. For a given interaction strength, 
bosonic bound states always have a larger velocity of spreading than their fermionic 
counterparts. When the interaction strength is increased, the spreading velocity 
decreases continuously from its free one-body value down to zero or to a much smaller 
limiting value. 

A natural extension of the present work consists in considering the quantum walk 
performed by bound states of more than two identical particles. A classical analogue 
consists of multi-pedal molecular devices whose legs perform random walks, known as 
molecular spiders [47]. Their behavior has been investigated theoretically, both on a 
one-dimensional [48] and a two-dimensional [49] substrate. The quantum version of 
the problem yields a special kind of A^-fermion bound state, which can be analyzed 
by techniques from integrable systems. This will be the subject of future work [50]. 
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Appendix A. Stationary properties of a quantum system 

In this appendix we investigate the stationary properties of a finite quantum system 
from a very general standpoint. Consider a quantum system whose Hilbert space has 
hnite dimension A, working in a preferred basis ja) (a = 1,...,A). In this basis, 
the Hamiltonian is given by an A x A Hermitian matrix H. Assume that the energy 
eigenvalues (n = 1,..., A) are non-degenerate. Let jn) be normalized eigenvectors, 
so that H\n) = A„jn). Assuming the system is initially in state Ja), we have 

IV'W) = ^e“‘^"*ln)(nla). 


n 


(A.l) 



Interacting quantum walkers: Two-body bosonic and fermionic bound states 


29 


The probability of observing the system is state \b) at time t reads therefore in full 
generality 

Pab{t) = |(6|V’(0)P = ^e‘(^”“-®"‘^‘(6|m)(TO|a)(a|n)(n|6). (A.2) 

m,n 


We are interested in the stationary transition probabilities 

1 f* 

Qab = lim - / Pab{t')dt'. 
t^OO t Jq 

In the absence of spectral degeneracies, only diagonal terms (m 
sum in (A.2) contribute. We thus obtain 

Qab = ^|(a|n)n(&|n)p. 


(A.3) 

n) in the double 
(A.4) 


The matrix Q is real symmetric and positive definite. It reads indeed || 
Q = RR^, 


(A.5) 


with 


Ran = \{a\n)\‘^. (A.6) 

The matrix Q is non-trivial in general. This is a manifestation of the well-known 
fact that a finite isolated quantum system does not equilibrate, in the sense that its 
stationary properties remember its initial state forever. 

We now turn to the explicit example of a confined quantum walker, i.e., a tight- 
binding particle on a finite segment of N sites labelled a = 1,... ,iV. The preferred 
basis is chosen to be local in space. The Hamiltonian reads 


(a|i?|'i/') = V'a-l-l + V'a-1, 


(A.7) 


where (a|^/’) = ipa and with Dirichlet boundary conditions -tpo = f’N+i = 0. We have 


En = 2 cos 


rnr 

N + V 


{a\n) = 


2 . amr 

■ sm ■ 


N+1 N+V 


(n = 1,..., N), and so 
Qab — 


E 


sm 


bmr 


N + 1 N+1 


■ sm 


^N+1 

This sum can be worked out explicitly. We thus obtain 

I 


Qab — 


N+1 


1 

1 + 2^a,b + 2^a-i-b,N+l 


(A.8) 


(A.9) 


(A.IO) 


With respect to its uniform background value I/(A^ + 1), the stationary 
probability Qab is thus enhanced by a factor 3/2 both at the starting point (6 = a) 
and at the symmetric position {b = N + 1 — a). In the particular situation where 
the starting point is the middle of an odd segment {N odd and a = {N + l)/2), the 
enhancement factor of the return probability reaches 2. 

The stationary mean value of the position A of a walker launched at site a, i.e., 

{X)=Y,bQab=^^, (A.ll) 

b 


denotes the transpose of the matrix R. 
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is dictated by symmetry, and therefore independent of the initial state. The 
corresponding variance, 

{X^)-(Xf =J2b^Q^,- 

b 



N{N-1) , (iV + 1 - 2 o) 2 
12 4(fV+ 1) 


(A.12) 


however shows a dependence on the initial position a of the walker. 
For iV = 4 we obtain 



/3 2 2 3\ 
2 3 3 2 1 
2 3 3 2 
V 3 2 2 3/ 


(A.13) 


The first row, i.e., 

3 1 

Qii = Qi4 = -pr, Qi2 = Qid, = 77 (A.14) 

10 5 

agrees with the result (3.34) of a full dynamical analysis of the two-fermion bound 
state with N = £ = 4,. 
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